0(N) symmetry-breaking quantum quench: Topological defects versus quasiparticles 
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We present an analytical derivation of the winding number counting topological defects created by 
an O(N) symmetry-breaking quantum quench in N spatial dimensions. Our approach is universal in 
the sense that we do not employ any approximations apart from the large- AT limit. The final result 
is nonperturbative in N, i.e., it cannot be obtained by an expansion in 1/N, and we obtain far less 
topological defects than quasiparticle excitations, in sharp distinction to previous, low-dimensional 
investigations. 



I. INTRODUCTION 

In contrast to the vast amount of literature regarding 
static properties (e.g., universal scaling laws) of phase 
transitions - both thermal and at zero temperature - 
we are just starting to understand their dynamical fea- 
tures, especially the behaviour during a time-dependent 
sweep (quench) through the critical point. This topic has 
attracted increasing interest in recent years, see, e.g., [l|- 
[9| and, again, some universal properties became evident 
[10(. For example, during a symmetry- breaking (second- 
order) dynamical phase transition, the diverging response 
time inevitably entails nonequilibrium processes and so 
the initial quantum (and thermal) fluctuations are am- 
plified strongly, ultimately determining the final order 
parameter distribution. If the final phase permits topo- 
logical defects (e.g., vortices in superfluids) , they will 
generally be created in such a quench via (the quantum 
version of) the Kibble-Zurek mechanism. The latter oc- 
curs in many diverse physical settings, for instance, in 
nonequilibrium phase transitions during the early uni- 
verse (lTl - [l3| or in condensed matter systems fl4| . 

Unfortunately, due to the inherent computational com- 
plexity of such scenarios, explicit calculations are difficult 
in general, and thus often rather uncontrolled assump- 
tions and approximations (e.g., Gaussianity [til. fl5|) have 
been invoked. For example, the correlation function after 
the transition has been used to infer the number of cre- 
ated quasiparticle excitations (see, e.g., [III). The quasi- 
particle number is, then, supposed to directly yield an 
estimate for the topological defect densities generated 
by the quench. For special cases such as the (exactly 
solvable) one-dimensional quantum Ising model (where 
the only excitations are topological defects, i.e., kinks 
P, [U, I!?]]), such an approach might give the correct an- 
swer - but in general, this will not be the case, as we will 
argue below. 

In the following, we consider a rather general O(N)- 
symmetry breaking quantum quench and study the cre- 
ation of topological defects (hedgehogs in the case consid- 
ered) via calculating their winding number. In order to 
base our derivation on a well-defined expansion, we con- 
sider the large- N limit. Apart from the large- A^ limit, no 



further approximations will be needed, i.e., our results 
will be quite universal. Moreover, similar to analogous 
large- A^ approaches in condensed matter and field theory 
(assuming that there is no critical value of N where the 
system changes drastically) , we expect our results to ap- 
ply qualitatively also to finite N (e.g., N = 3), which are 
accessible to experimental tests. Bose-Einstein conden- 
sates, in particular, permit the time-resolved observation 
of the defect formation mechanism due to the compar- 
atively long req-equilibration time scales of these dilute 
ultracold quantum gases [lj| . 



II. EFFECTIVE ACTION 

As a first step, we construct a general effective action 
for an 0(A^)-model in terms of the A^-component field 
<fi — (0i, • ■ • , 4>n), which determines the order parameter. 
To this end, we start from the equation of motion with 
an arbitrary function / 



</> = /(0, 0, v 2 0, v 2 ^, v 4 0, v 4 0, . . . ) 



(i) 



In order to avoid run-away solutions and to facilitate a 
proper quantum description, we have assumed the ab- 
sence of time derivatives of third or higher order. The 
initial state (before the transition) obeys the O(N) sym- 
metry: {(p a ) = and {4> a {x)4>b{x')) oc 5 ao , etc. As stated, 
in all of our calculations, we employ the large-A^ limit 
assuming N ^> 1. In this case, O(N) invariant combi- 

nations such as </> = 4>\ + • • ■ + tp 2 N are sums of many 
independent quantities on an equal footing [20| . Consid- 
ering commutators of such combinations, we obtain the 
well-known fact that their leading contribution (in the 
large- A^ limit) behaves as a c-number whereas the (clas- 
sical and quantum) fluctuations scale with \/N (cf. the 
law of large numbers). Therefore, we may approximate 



4> 2 = (4> 2 )+o(Vn), (4> 2 }=o(n), 



(2) 



arriving at a semi-classical (mean-field) expansion valid 
in the large- Af limit. As a result, we may approximate 
the nonlinear terms in the equation of motion ([lj, for 
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~3 -2 - 

example by <f> ss (0 )<p, arriving at a linearized descrip- 
tion. This leads us to the most general linear and local 
O(N) invariant effective action containing up to first time 
derivatives of the fields <fi = (01: ■ • ■ , 4>n) 



(3) 



with arbitrary Fourier space functions F(k 2 ) and G(k 2 ). 

III. PHASE TRANSITION 

From Eq. ©, we derive a Klein-Gordon type disper- 
sion relation [to 0(k 2 )] for the linearized fluctuations, 



0J 2 {k) 



G(k 2 
F(k 2 



2 4 

= m c 



c 2 k 2 



0(k 4 



(4) 



Initially, all modes are stable, u> 2 (k) > 0, since we lin- 
earize around the initial [0(A r )-symmetric] state. Af- 
ter the 0(iV)-symmetry breaking transition, however, 
the state (0) = is no longer stable and the system 
"wants" to roll down to a state with (</>) ^ 0. Typ- 
ically (for second-order transitions, i.e., without meta- 
stability), this implies that some of the modes become 
unstable, w 2 (fc) < 0, cf. Fig.[TJ Since Eq. ^ is already a 
result of the large- N limit, we assume that uj 2 (k) is inde- 
pendent of N 3> 1 (otherwise the group and phase veloc- 
ities would either diverge or vanish in the limit N — > oo). 
Furthermore, modes with sufficiently large k should be 
stable Lu 2 (k f oo) > 0, so that the unstable interval in 
which ui 2 (k) < is assumed to be finite. 

So far, our results were independent of the number D 
of spatial dimensions due to isotropy. In the following, we 
set N = D in order to facilitate the creation of topological 
defects in the form of hedgehogs (see below). So strictly 
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FIG. 1. Two generic examples for the evolution of the disper- 
sion relation J3| during a symmetry-breaking phase transi- 
tion, see also [101 ] . Initially (dotted line), all k- values are sta- 
ble, uj 2 (k) > 0. At the critical point (dashed line), the disper- 
sion relation touches the fc-axis, and after the transition (solid 
line), modes in a finite fc-interval become unstable, u> 2 (k) < 0. 
The left panel corresponds to a case where u> 2 (k = 0) = m 2 
in Eq. Q remains positive while c 2 changes sign (see, e.g., 
[2lT|h whereas, in the right panel, u) 2 (k = 0) = m 2 becomes 
negative. In both cases, however, there is a dominant wave 
vector fc* (for large N), as indicated by the vertical arrows. 



speaking, we consider the simultaneous limit N — > oo and 
D — > oo and assume that these limits commute. Since 
all relevant quantities such as ((^air, t)(f> D (r', t)) depend 
on \r — r'\ only (isotropy), the large- £> limit basically 
just affects the integration measure d D k, which strongly 
supports this assumption. With a Fourier expansion of 
Eq. ([5]) , we obtain the two-point function after the quench 



Sab 



(2n) N / 2 



(- r >) 



with L = \r — r'\. The Bessel functions J„ with in- 
dex v = N/2 — 1 arise from the integration over all fc- 
directions and the factors Cjjr and depend on the ini- 
tial state (for example the temperature) as well as quench 
dynamics, and are roughly independent of N . As ex- 
pected, we obtain an exponential growth of the unstable 
modes, which have uj\ < 0, after the phase transition - 
which then seeds the creation of topological defects. Of 
course, due to the growing modes, the linearization in 
Eq. j3]) will fail eventually - but for N f oo, the time 
t until which the linearization and thus Eq. ([S]) applies 
does also grow. Therefore, we may distinguish basically 
three phases following the quench: First, we get a period 
of exponential growth of the modes where the linearized 
description in Eqs. ^ till ([5]) applies. Then, nonlin- 
ear effects set in and lead to a saturation of this growth 
and possibly an oscillation around the new energy min- 
imum. Finally, the topological defects created by the 
quench start to "feel" the attraction between hedgehogs 
and anti-hedgehogs leading to their approaching each 
other and eventual annihilation. This general picture 
has been qualitatively confirmed by numerical simula- 
tions for N = 2. The topological defects are seeded 
in the first phase (exponential growth) and slowly disap- 
pear in the final phase. Therefore, we expect to obtain 
a good estimate for the maximum number of created de- 
fects at long length scales and intermediate times from 
our linearized analysis. 

In addition to the exponentially growing (in t) modes 
at finite k, the integral ([S]) does also yield a huge contri- 
bution from large k, because the phase space factor fc^ -1 
rapidly rises with k for large N. This gives rise to a strong 
UV singularity of the two-point function oc |r — r'\~°( N \ 
see the discussion below of regularizing this UV diver- 
gence. 



IV. TOPOLOGICAL DEFECTS 

After the symmetry-breaking transition, the ground 
state is degenerate and can be specified by a nonvanishing 
expectation value (cf)) ^ 0, which singles out a preferred 
direction given by the unit vector n = (4>)/\ (cf>) |. Thus 
the original O(N) symmetry is broken down to 0(N — 1), 
i.e., rotations around the n-axis, and the ground-state 
manifold corresponds to the surface Sn-i of a sphere in 
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TV dimensions 0(N)/0(N — 1) ~ <Sjv-i- Remembering 
the homotopy group 7Tjv-i (<S/v-i) = Z, we see that topo- 
logical point defects in the form of hedgehogs exist in 
N spatial dimensions (22l . |23| . These defects correspond 
to nontrivial mappings from the ground-state manifold 
Sn-i onto the surface Sn-i of a sphere in real space, 
characterized by a winding number 9t € Z, which reads 

0,m 



m = 



£abc...£ 



a/37-. 



r(iv) ||<Sjv-i|| 



dS a n a (d p n b )(d^n c )..., (6) 



where || «Sjv-i lh 2n N / 2 /T(N/2) is the surface area of the 
unit sphere in N dimensions. Starting with the O(N)- 
symmetric state as the initial state, we cannot simply 
insert n = (<p)/\(<p)\ since (<p) vanishes. Therefore, we 
use a quantum operator n instead, which must be defined 
appropriately, and allows for a derivation of the probabil- 
ity distribution of the quantum winding number 91 in a 
given volume from the above general expression. In par- 
ticular, the expectation value of the winding number is of 
course zero, (OT) = 0, but its variance (9T 2 ) is in general 
not. Setting directly h oc </>, we see that the variance 
(9r) is plagued with UV divergences similar to other 
quantities containing products of quantum fields at the 
same space-time point. This is due to the fact that quan- 
tum fluctuations of cj> at arbitrary fc-scales [cf. the strong 
UV singularity mentioned after Eq. ^\ would in general 
contribute to (9r). In fact, even the 0(AT)-symmetric 
initial ground state can be viewed as a "quantum soup" 
of virtual hedgehog-anti-hedgehog pairs which are con- 
stantly popping in and out of existence. Here, we are 
not interested in those virtual short-lived defects, but in 
long-lived hedgehogs, which are created by the quantum 
quench. Therefore, we have to insert a time-averaged 
unit vector defined via 



n(r) 



1 

~Z 



dtg(t)<j>(t,r), 



(7) 



with a smooth smearing function g(t) and the normal- 
ization Z = ([J dtg^^t.r)] 2 ) 1 / 2 + 0(VN), where we 
have used the large- iV (mean-field) expansion. This time- 
average now suppresses all (rapidly) oscillating modes 
with w 2 , > and only leaves the growing modes u> 2 < 0. 
After this UV-regularization, the integral in the two- 
point function (|S|) will be dominated by only a few modes 
in the vicinity of a certain wavenumber fc»: In view of 
the phasespace factor fc^ -1 in the dominant contri- 
bution for large N [3(| will arise from the largest k value 
for which uj\ < 0, i.e., close to the zero of w 2 ., cf. Fig.Q] 
Thus, we can evaluate ([5]) in saddle point approximation 
and obtain the correlator for the time-averaged direction 
vector in ([7]) 

(n a (r)n b (r')) = 2" ^+ 1) J -0^ 6 ab = f(L)S ab ,(8) 

where higher-order terms cx iV~ 3 / 2 stemming from the 
normalization Z in ((7]) have been omitted. These higher- 



order terms vanish for N — s- 00 so that ((SJ becomes in- 
deed exact and, in view of the asymptotic behavior of the 
Bessel functions [22|, a Gaussian correlator 



ffn 1 / k 2 L 2 



(9) 



follows. Thus, the typical linear domain size (correla- 
tion length) is given by L C orr = 0(y/N/k*). At ex- 
tremely large distances L = 0(N/k*) (where the first 
nontrivial zero of the Bessel function J v is located) , there 
are oscillatory deviations, but in this regime, the cor- 
relator is already exponentially small. We emphasize 
that the Gaussian form of f(L) stems from the large 
TV limit of the exact expression in ([8]), and is not as- 
sumed a priori. Furthermore, as may already be ob- 
served in Eq. ([SJ , the emergence of a dominant wavevec- 
tor &;* implies the cancellation of all time-dependence, 
i.e., the time-dependence of the growing part of ap- 
proximately separates such that the time-averaged unit 
vector (J7J becomes independent of g(t) and thus station- 
ary (in the regime under consideration). So the emer- 
gence of a dominant scale in the correlator, fc*, a poste- 
riori justifies the introduction of the UV regulator g(t), 
which only affects the rapidly oscillating modes at larger 
k but not the observables we are interested in. 



V. SCALING LAWS 

Now we are in a position to derive the dependence of 
(9r) on N and the enclosed volume. Inserting Eq. (J7J 
into the winding number variance (9r) from Eq. ([5]), 
we obtain the expectation value of the product of 2N 
fields (j> a , which factorizes into N two-point functions (JSJ. 
Since these functions are completely regular, we may ap- 
ply Gauss' law to the two surface integrals occurring in 
(9l 2 ), and get after some algebra (28| 



m z = 



NNl 
|«5jv-i| 



d N rd N i 



1 



d 



L N_1 dL 



2i 

dL 



N 



(10) 



For a sphere of radius R, V — {r : r 2 < R 2 }, we can 
evaluate this expression and finally obtain a single inte- 
gral of the form 



m z ) = — R 

7T 



it/2 

N ' d9 



N 



cos9^-(2Rsm9) 







Note that these general expressions (fTU|) and pT|) for (9r) 
are neither restricted to the Bessel functions ((SJ), nor to 
large N, but might be employed for any correlator of the 
form (h a (r)hb(r)) — 5 ab f(L) in any dimension N > 2. 
Let us discuss the scaling of (CH 2 ) with respect to R and 
N 3> 1, using the Gaussian correlator ©. For radii far 
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above the correlation length R L corr = viV/fe*, the 
winding number variance behaves as 



= 



-3/2 



KR 



1 + 0(1/VN) 



N-l 



(12) 



We observe that (9r) scales with the area R N ~ X of the 
hyper-surface enclosing the defects. Apart from the pref- 
actor e _3 / 2 /v A /V, this area scaling is quite universial as 
it holds for spherical volumes in any dimension N > 2 
- provided we assume short-range correlations - and can 
already be inferred from Eq. ([5]): If we calculate (*Tt 2 ) 
using ([5]), we obtain two hyper-surface integrals. Due to 
isotropy. the first one yields R N_1 while the second inte- 
gral averages over the distance \r — r'\ between the two 
points on the surface. Assuming short-range correlations 
only, this second integral becomes independent of R (for 
large R) and gives h{N)k^~ x with some function h(N). 
Note, however, that the assumption of short-range cor- 
relations is crucial and nontrivial in this argument: For 
vortices in two dimensions, for example, we obtained log- 
arithmic corrections to the "area" scaling, (9r) oc RlnR 
[H, since the correlator fell off quite slowly at large L. 

If the radius R shrinks and approaches the correlation 
length R ~ L corr = 0{-\/~N the winding number 
variance decreases rapidly (for N 3> 1). A sphere with 
R = 0(\f~N /k*) would then contain around one defect 
(or anti-defect) on average, (9r) = 0(1), which deter- 
mines the total defect density. For even smaller radii, far 
below the correlation length R <C Lcorr = 0(y/N/k*), 
the above formulae would yield a scaling (*Jl 2 ) ~ R 2N , 
i.e., an exponential suppression (for large N). However, 
the precise functional form (fTTj) should not be trusted 
upon in this regime since we have neglected 0(l/y/~N)- 
corrections in our derivation, which is problematic if the 
final result is exponentially small. From a more physical 
point of view, the mean-field approximation @ breaks 
down near the core of a defect (where n becomes ill- 
defined), which renders Eq. (fTTj) questionable for too 
small volumina. For small R, one would expect a volume- 
type scaling (9r) ~ R , i.e., the typical behaviour for 
uncorrelated defects, which should be the case if there 
is one hedgehog at most. The exponential suppression 
(m 2 ) ~ exp{-0(iV)} for large N and small R < L corr 
should still be correct, as this just reflects the diminishing 
probablity of reversing field orientation in all directions 
when increasing N. 

The area scaling (fT2"j) of the net defect number (*Xt 2 ) 
can be interpreted as the occurrence of a confined phase 
of bound defect-antidefect pairs. Only pairs where one of 
the partners is contained within while the other is out- 
side the integration volume would yield net winding num- 
ber, whereas those pairs entirely inside or outside do not 
contribute. Hence a scaling with surface area instead 
of volume is natural for short-ranged correlations. On 
the other hand, there could, in principle, also exist a 
de-confined phase of quasi-free hedgehogs similar to the 
quark-gluon plasma of quantum chromodynamics. Such 



a de-confined phase might occur if defect density and 
temperature are sufficiently high and any bound pairs are 
broken up again by thermal quasi-particles. In that case, 
defects and antidefects would be randomly distributed 
and a volume scaling of the winding number variance 
(01 2 ) follows. Since we did not make any assumption in 
our derivation apart from the large N limit (where the 
results become exact), we can clearly distinguish between 
the two phases (confined or de-confined). 



VI. STATISTICS 

In a similar manner, we can calculate the higher mo- 
ments of the winding number. Again, by exploiting the 
fact that we have short-range correlations, the large- R 
limit of the next nontrivial moment can be inferred from 
pure combinatorics in the analysis of the four integrals 
occurring in 

<0T 4 ) = 3(<Tt 2 ) 2 + 0(R N - 1 ) = 0(R 2N - 2 ) . (13) 

Analogously, the leading terms of ( < Jt 2n ) are given by 
(2n-l)!!(^l 2 )" with (2n-l)!! = (2n-l)(2n-3) . . .5-3. 
For large R, the winding number CI £ Z can be ap- 
proximated by a continuous variable OT g R and thus 
its full statistics is given by the inverse Mellin transform 
of (2n- 1)!! = 2T(n+l/2)/V7r, which yields the Gaus- 
sian probability distribution p(9T) cx cxp{— 7 2 9t 2 }, with 
I/7 2 = 2(m 2 ). We note that, like in Eq. ©, the Gaus- 
sianity is not assumed but derived from first principles in 
a given limit - for small R (small 91), for example, there 
will be deviations from a Gaussian distribution. 



VII. CONCLUSIONS 

Based on a very general 0(A^)-invariant effective ac- 
tion, we presented an analytical derivation of the wind- 
ing number counting the defects created by a symmetry- 
breaking quantum quench in the large- N limit. Consis- 
tent with previous calculations jUJ, our result (fT!?)) is 
nonperturbative, i.e., it does not admit a Taylor expan- 
sion in 1/N. As another result, we find that the typi- 
cal distance between defects scales with the correlation 
length 0{sfN /fc„). By contrast, the typical distance be- 
tween quasiparticle excitations (e.g., Goldstone modes) 
does not increase with N. This can be understood by re- 
calling that the total energy of the system (which scales 
with N) in a given volume has to be distributed among all 
the quasiparticle excitations, whose typical energy is de- 
termined by the dispersion relation ui 2 (k) and thus inde- 
pendent of N. Therefore, we conclude that the quasipar- 
ticle spectrum alone does not yield any direct information 
about the generation of topological defects in general. 
This situation is quite different in the one-dimensional 
quantum Ising model, where topological defects (kinks) 
are the only quasiparticle excitations [l], which fre- 
quently led to the assumption in the literature that this 
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is generic. We demonstrated here that identifying quasi- 
particle excitation and defect numbers created by a quan- 
tum quench can be quite misleading. 

The crucial difference between quasiparticles (whose 
number can be derived via a perturbative expansion in 
1 /N) and topological defects (which are nonperturbative) 
can be illustrated by the following intuitive picture: Con- 
sidering a discrete regular lattice with a unit direction 
vector rii at each lattice site i, a quasiparticle excitation 
occurs if rii ^ rij for two neighbours A topologi- 
cal defect at the site i, one the other hand, means that 
the unit vectors rij of all neighbouring sites either point 
away or towards the site i. For large N, this is obviously 
a much stronger condition. 

It is also worth noting that the derived area scaling 
(9r) oc R 1 ^^ 1 is inconsistent with the random defect 
gas model (where defects and anti-defects are distributed 
randomly in the sample volume, corresponding to a de- 
confined phase) since this model would predict a volume 
law, i.e., i? w -scahng. We remark in this connection that 
the area scaling [26| (corresponding to a confined phase) 
we obtain can be interpreted by a random n-field model 
on the hyper-surface with the correlator Eq. ([5]). repre- 
senting a generalization of the random phase walk model 
for N = 2 (cf. the result of in which reasonable agree- 
ment with the experiment reported in fl8| was obtained) . 

Finally, we would like to stress that our result is quite 
universal, i.e., it is valid for very general dispersion re- 
lations of the O(N) model (cf. Fig.[I| and just relies on 



the large-iV limit without any further approximations. 
Moreover, as indicated below Eq. (fU2"j) . we expect that the 
general picture does still apply qualitatively for smaller, 
and thus experimentally accessible values of N, for ex- 
ample N = 3. In particular, this should be true for fast 
quenches, where we have a well-defined period of expo- 
nential growth of the unstable linear modes, while non- 
linear effects (saturation of this growth, oscillations, and 
finally defect annihilation, see [23) occur much later. In 
this case, one may find (instead of 1 /N) another small pa- 
rameter (e.g., the diluteness of the gas) in order to moti- 
vate the underlying effective action in analogy to Eq. ^ . 
For N ^> 1, universality will be partially lost and the 
dependence on the dispersion relation, for example, will 
be stronger. For instance, it might then be necessary to 
introduce a time-dependent critical fc* = fc*(t), which is 
not close to the zero of w 2 (fc), but near the actual mini- 
mum of oj 2 (k). 
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